Pattern Formation as a Signature of Quantum Degeneracy in a Cold Exciton System 
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The development of a Turing instability to a spatially modulated state in a photoexcited electron- 
hole system is proposed as a novel signature of exciton Bose statistics. We show that such an 
instability, which is driven by kinetics of exciton formation, can result from stimulated processes 
that build up near quantum degeneracy. In the spatially uniform 2d electron-hole system, the 
instability leads to a triangular lattice pattern while, at an electron-hole interface, a periodic Id 
pattern develops. We analyze the mechanism of wavelength selection, and show that the transition 
is abrupt (type I) for the uniform 2d system, and continuous (type II) for the electron-hole interface. 
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PACS numbers: 

Recently, striking spatial photoluminescence (PL) pat- 
terns Q, S S Hi which span macroscopic scales in 
excess of 100 /im, have been observed in photoexcited 
AlGaAs/GaAs quantum well (QW) structures. In addi- 
tion to concentric rings and 'bright spots', the electron- 
hole system exhibits an abrupt transition at ca. 2 K 
in which the outermost ring 'fragments' into regularly 
spaced beads of high PL intensity 0,0] . While the gross 
features of PL have been explained within classical frame- 
work, attributing the internal rings to nonradiative exci- 
ton transport and cooling and the outermost rings and 
'bright spots' to macroscopic charge separation 0,0, 0|, 
the origin of the instability remains unidentified. 

Spatially modulated exciton density is not to be ex- 
pected in QW system designed so that excitons interact 
repulsively as electric dipoles and thus do not form 
droplets |6(. The macroscopic character of ring frag- 
ments, of 10 — 30 /im in size and containing about 10 
excitons each, and the abrupt temperature dependence, 
call for an explanation involving a symmetry-breaking 
instability of a homogeneous state to a patterned state. 
Such behavior is reminiscent of the instability predicted 
by Alan Turing to occur in a reaction-diffusion sys- 
tem. The Turing instability is known to occur in certain 
chemical reactions 0, 0, , and is also believed to be 
relevant for pattern formation in biological systems ^lj • 

In this work we propose a mechanism, based on the 
kinetics of exciton formation from optically excited elec- 
trons and holes, that can lead to instability in the exciton 
system. Interestingly, these kinetic effects become espe- 
cially strong in the regime near exciton quantum degener- 
acy, due to stimulated enhancement of the electron-hole 
binding rate. The transition to a state with a spatially 
modulated exciton density reveals itself in the spatial PL 
pattern, and presents a directly observable signature of 
degeneracy. 

Although, in itself, an observation of an instability does 
not constitute unambiguous evidence for degeneracy, it 
may complement other manifestations discussed in the 



literature, such as changes in exciton recombination [l2T | 
and scattering 0] rates, in the PL spectrum Q, absorp- 
tion 0| and PL angular distribution 0]. While linking 
the observed instability with degeneracy is premature, 
our main aim here is to present the Turing instability 
from a broader viewpoint, as a novel and quite general 
effect of quantum kinetics that can help to identify the 
regime of Bose-Einstein condensation (BEC). 

Here we consider a transport theory 0] formulated in 
terms of electron, hole, and exciton densities n ej h, x obey- 
ing a system of coupled nonlinear diffusion equations: 
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The nonlinear couplings account for exciton formation 
from free electron and hole binding. Here, no attempt 
has been made to describe the detailed and complicated 
density and temperature dependence of the physical pa- 
rameters entering the model, nor to account for nonequi- 
librium exciton ener gy distribution and cooling due to 
phonon emission |18Lll7l |. Instead, we adopt a more phe- 
nomenological approach and assume that, as a result of 
cooling, the system can be described by an effective tem- 
perature, which leaves the densities n e; h, x as the only im- 
portant hydrodynamical variables. The sources J e and Jh 
in Eq. (JJJ describe the carrier photo-production, as well 
as the leakage current in the QW structure. 

In general, one can expect the electron-hole binding 
rate w and, to a lesser extent, the exciton recombination 
rate 7 to depend sensitively on the local exciton density 
n x . Of the several mechanisms that could lead to such a 
dependence at low temperatures close to exciton degen- 
eracy, perhaps the most important in the present context 
involve stimulated electron-hole binding processes medi- 
ated by phonons. These processes enhance the binding 
rate by a factor / = 1 + N^, where denotes the oc- 
cupation of exciton states. In thermal equilibrium, and 
at low temperatures, one can ignore the reverse processes 
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of exciton dissociation: The binding energy, carried away 
by phonons, is much larger than k^T . 

For a degenerate exciton gas with Ne=q > 1, the dom- 
inant process involves scattering into the ground state 
and the stimulated enhancement factor is expressed as 

2gm x k^T 
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Here to x ~ 0.21 mo represents the exciton mass, and 
g denotes the degeneracy (in the indirect exciton sys- 
tem, the exchange interaction is extremely weak, and 
g = 4). Equivalently, when reparameterized through 
its dependence on temperature, u = Tq/T where Tq = 
(irH 2 /2<?to x /cb) n x is the degeneracy temperature. At 
T ~ Tq (equivalently n x ~ rio), there is a crossover 
from classical to quantum Bose-Einstein statistics, and 
the stimulated enhancement factor / increases sharply. 

Qualitatively, the stimulated transition mechanism for 
hydrodynamic instability can be understood as follows: 
A local fluctuation in the exciton density leads to an in- 
crease in the stimulated electron-hole binding rate. The 
associated depletion of the local carrier concentration 
causes neighboring carriers to stream towards the point 
of fluctuation presenting a mechanism of positive feed- 
back. The wavelength, determined by the most unstable 
harmonic of the density, characterizes the lengthscale of 
spatial modulation in the nonuniform state. 

Before turning to the analysis of instability, it is use- 
ful to discuss intrinsic constraints on the dynamics |JT|| 
due to electric charge and particle number conservation. 
These are obtained by considering the linear combina- 
tions (e) — (h), (e) + (h) + 2(x) of the transport equa- 
tions In both cases, the nonlinear term drops out 
and one obtains linear equations 



L e n e — Lh^h = Je — Jh 

L c n c + L h n h + 2L x n x + 2^n x = J c 
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Note that, since the ori- 



gin of the relations il-il II) is routed in conservation laws, 
they are robust and insensitive to the exact form of the 
electron-hole binding term. 

Initially, let us consider a system in which the sources 
are constant J e (r) = Jh(r) = J, realized by a spatially 
extended photoexcitation. In this case, ignoring the de- 
pendence of the recombination rate 7 on density, we have 



= J/j, n e , h = (J/w{n x )) 1/2 . 



(5) 



The stability of the system can be assessed by linearizing 
Eqs. Q about the uniform solution JSJl with a harmonic 
modulation <5n Gj h, x °c e e . Using (l.'il 111 and writing 
L e Sn e = LhSrih — —L x Sn x — "fSn x , one can express 5rih,x 
in terms of 5n c and obtain 



L e (A,k)+ 7 ^ 
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where u = d In w/d In n x is evaluated at the steady state 
©, and L c(hiX) (A,k) = A + D c(hyX) k 2 . Solving Eq. ©, 
one obtains the growth rate dispersion A(k) from which 
one can infer from the stability criterion, Re A < 0, that 
the system becomes unstable when 
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where l x = W D x /^ denotes the exciton diffusion length 
and r = Z? X (ZA7 1 + )n x /n e . Eq. (7J has solutions if 
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as illustrated in Fig. ^ At u — u c , we obtain the most 
unstable wavenumber fc* — r 1 / 4 / ' l x selected by competi- 
tion of the stimulated binding and diffusion processes. 

The binding rate w(n x ) builds up near degeneracy due 
to the growth of stimulated processes, leading to insta- 
bility at temperatures approaching Tbec- (For a weakly 
interacting Bose gas the transport coefficients, and thus 
the constant r, are practically insensitive to the degree 
of exciton degeneracy at T > Tbec-) 
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FIG. 1: Graphical solution of Eq. Q that selects the most 
unstable wavelength. Inset: The 3-fold symmetric star of 
wavevectors describing the modulation near the instability 
threshold and the corresponding triangular pattern of exciton 
density variation. 

To what extent are these results insensitive to the ori- 
gin of the nonlinearity in the binding rate? If enhanced 
by intraband Auger processes, which transfer the binding 
energy released in exciton formation to other excitons, 
one expects the binding rate w to scale linearly with lo- 
cal exciton density, viz. w(n x ) = u>o(l+n x /no), where no 
denotes some constant involving a ratio of the two-body 
and three-body cross-section of the electron and hole in 
the presence of excitons. Crucially, in this case, the left 
hand side of Eq. (jHJ is bounded by unity, while the right 
hand side is in excess. Therefore, at least over the param- 
eter range considered here, one can infer that a simple 
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linear scaling of the binding rate with density does not 
lead to instability. Indeed, the instability may be used to 
discriminate against certain mechanisms in the kinetics 
of exciton formation. 

Turning to the discussion of the spatial pattern result- 
ing from the instability, we note that the wavevector se- 
lection determines its modulus, but not direction. At 
threshold u — u c all modes with |k| = fc* become un- 
stable simultaneously. The resulting 2d density distri- 
bution can be found by considering the effect of mixing 
different harmonics due to higher order terms in ex- 
panded in <5n e , x about the uniform state. Since these 
equations contain quadratic terms, the favored combi- 
nation of harmonics is a 3-fold symmetric star k^' = 
k*(cos( 2 fj+6),s'm(^j+8)), j = 1, 2, 3, with the param- 
eter 6 describing the degeneracy with respect to 2d rota- 
tions. This leads to a density distribution Sn cx J2j e ±kj ' r 
with maxima arranged in a triangular lattice. 

On symmetry grounds, since the triangular lattice pat- 
tern is stabilized by quadratic terms, the mean field anal- 
ysis predicts that the transition to the modulated state in 
this case is abrupt, of a type I kind. Indeed, the triangu- 
lar lattice geometry arises in various 2d pat tern selection 
problems, from Benard convection cells [lj| to the mixed 
state of type II superconductors 

To explore the application of these ideas to the Id mod- 
ulation seen in exciton rings Q, Q , one must first deter- 
mine the profile of the uniform distribution. The rings 
represent an interface between regions populated by elec- 
trons and holes at which they bind to form excitons. The 
steady state is maintained by a constant flux of carriers 
arriving at the interface. The parameter regime which is 
both relevant and simple to analyze is that of long ex- 
citon lifetime 7 -1 where the diffusion length i x exceeds 
the range of the electron and hole profile overlap. In this 
case, approximating the source of excitons by a straight 
line c5(x), where c is the total carrier flux and x is the 
coordinate normal to the interface, the exciton density 
profile is given by (c£ x /2D x )e~' x '' ix . Accordingly, one 
can seek the electron and hole profile treating w(n x ) as 
constant and restoring its dependence on n x later when 
turning to the instability. The profiles can be inferred 
from two coupled nonlinear diffusion equations 

D c ( h )dln c(h ) = wn c n h , (9) 

with the boundary condition: Dcd x n^\± oc = ±c9(±x). 
From Eq. @ one obtains D c n c — D^n^ — cx which allows 
the elimination of n^. Applying the rescaling n e (h) = 
C ^5e(h)/-C ) e(h)j where i = (DgD^/wc) 1 ^ 3 , one obtains 

dj g e = ffc(ffc - x), x = x/t. (10) 

From the rescaling one can infer that the electron and 
hole profiles overlap in a range of width I ~ c -1 / 3 while 

g e (\x\ » t) = x9(x) + O (|£r 1/4 e~ 2|5;|3/2/3 ) • 
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FIG. 2: Numerical solution of the (normalized) bound state 
wavefunction ip(x) representing fluctuation in the electron 
density, together with the effective potential U c g = g c + f]h — 
a o ffcffh of the Schrodinger-like equation ill II . 

Although Eqs. Q arc nonlinear, their diffusive charac- 
ter does not straightforwardly admit a spatial instability: 
A fluctuation in the position of the interface initiates 
an increased electron-hole flux which, in time, restores 
the uniform distribution. However, if one restores the 
dependence of the binding rate on exciton density, the 
same mechanism of positive feedback which character- 
ized the instability in the uniform system becomes ac- 
tive. To explore the instability, one may again expand 
linearly in fluctuations around the spatially uniform so- 
lution, g e (x, y) = g c (x)+8g c (x) e lky (similarly g h and g x ), 
where y is the coordinate along the interface and g e (x) 
denotes the uniform profile obtained from Eq. QlOf) . With 
£ x S> I, the exciton density remains roughly uniform over 
the electron-hole interface. Denoting this value by n x (0), 
in the vicinity of the interface, one may again develop 
the linear expansion w[n x ] ~ w[n x (0)] (1 + uSn x /n x (0)) 
where, as before, u = dlnw/d\nn x . Noting that Eq. @ 
enforces the steady state condition Sg e = Sg^, a lineari- 
sation of Eqs.Q) obtains the Schrodinger-like equation 

[-8j + {Ikf + g c + g h ] Sg c + r^rSsSh Sg x = 0, (11) 
together with the condition on the Fourier components, 

- S9M (l - Q^p) , (12) 

with Q 2 = (k£ x ) 2 + 1. Now, since the product g c g~h is 
strongly peaked around the interface, the typical con- 
tribution from the last term in (|llfl arises from Fourier 
elements ql ~ 1. Then, with £ x ^> £, the second contribu- 
tion to 5g x (q) can be treated as a small perturbation on 
the first and, to leading order, neglected, i.e. Sg x ~ — Sg c . 
In this approximation, the most unstable mode occurs at 
k = 0. Qualitatively, an increase in u will trigger an in- 
stability of the k = mode at a critical value u c when 
the linear equation first admits a non-zero solution for 
Sg c . At the critical point, the corresponding fluctuation 
in the electron density then acquires the profile of the 
(normalized) zero energy eigenstate ip{x). Numerically 
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one finds that the critical point for the instability oc- 
curs when u c /g x (0) — (2£/£ x )u c = a ~ 6.516, while the 
corresponding solution ip(x) is shown in Fig. El 

While the approximation above identifies an instabil- 
ity, a perturbative analysis of the ^-dependent correc- 
tions implied by l|12l) reveals that the most unstable mode 
is spatially modulated. To the leading order of perturba- 
tion theory, an estimate of the shift of u c obtains 
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dx g e (x)gh(x)ijj 2 (x) ~ 0.254^, and 



where a\ = J° 
1 f°° 

a>2{Q) = -f dxdx' g e (x)g h (x)ijj{x) 

^ J -oo 



With £ x » £, it will follow that Qt < 4, and the latter 
takes the constant value — 0.46D? 2 independent of k. 
Finally, minimizing 5u c with respect to fc, one finds that 
the instability occurs with a wavevector 

^1/3 



k c £ x ~ (a a 2 £ x /2£ i 



(13) 



implying a shift of u c by Su c /u c ~ (£/£ x ) 4 ^ 3 . As a re- 
sult, one can infer that the spatial modulation wavelength 
A c ~ i 1 / 3 ^ 3 is typically larger than the electron-hole 
overlap I, but smaller than £ x . Finally, an expansion 
of the nonlinear equations to higher order in fluctuations 
shows that, below the transition (i.e. for u > u c ), the am- 
plitude of the Fourier harmonic k c grows as (u — u c Y^ 2 . 
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FIG. 3: Amplitude of the spatial modulation of the exciton 
density g x along the electron-hole interface x — as a func- 
tion of the control parameter u/g x (0) as determined from the 
numerical solution of the nonlinear transport equation 1141 . 
A fit of the data to a square root dependence ~ (u — u c ) 1//2 is 
shown. Inset: Dimensionless electron g e and exciton g x den- 
sities for it/<; x (0) — 9.5 and £ x /£ = 4. Here periodic boundary 
conditions are imposed along the interface. 

Once the instability is strongly developed (or when 
£ x < I) the linear stability analysis above becomes unre- 
liable. Here one must turn to numerics. Having in mind 



the mechanism of stimulated scattering, Fig.|3|shows the 
results of a numerical analysis of the dimensionless non- 
linear steady-state equation 



Vx 9c = cxp 



5x(0) 



%c g c {g c - x), 



(14) 



where, using Eq. I|12|) . Sg x = g x — <7x depends non- locally 
on Sg c = g c — g c through the linear relation 



<5ffx(x) = -<^5c(x) + 
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and Ko denotes the modified Bessel function. Although 
(as in the experimental ring geometry) the modulation is 
constrained by the periodic boundary conditions imposed 
along the direction parallel to the interface, the critical 
wavenumber k c lies close to that predicted by Eq. ifH?)) . 
Similarly, the constraint leads to a value of u c /g x (0) a 
little in excess of that predicted by the linear stability 
analysis. Finally, the amplitude of the instability con- 
firms the square root dependence on (u — u c ) predicted 
by perturbation theory. 

In summary, we have shown that the realization of 
quantum degeneracy in a cold electron/hole-exciton sys- 
tem is signalled by the development of a spatial density 
modulation. Although our discussion is motivated by the 
photoexcited CQW system in which electrons and holes 
are spatially separated, the mechanism is quite generic 
applying also to geometries where the electron and hole 
sources J e ,h have a spatially independent profile. By con- 
trast, the instability mechanism appears to depend sen- 
sitively on there being a strongly nonlinear dependence 
of the electron-hole binding rate on the exciton density 
pointing to the importance of stimulated scattering. 
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